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ABSTRACT 


Submitted to the Department of Naval Architecture suid Marine Engineering 
on May 17» 1968 , in partial fulfillment of the requirements for the Mas¬ 
ter of Science Degree in Electrical Engineering and the Professional 
Degree, Naval Engineer. 

Previous electric-field-type electromechanical energy conveirters have 
relied on a variable capacitance, with terminal voltage constrained. An 
electrostatic generator is studied which employs a spatially varying charge 
distribution with the potential unconstrained. Energy is converted from 
mechanical to electrical form by means of a synchronous interaction between 
the excitation charge wave on the moving medium and the fixed load. Power 
output characteristics for both a continuum and discrete loading arrange¬ 
ments are derived and analyzed. A study of the effects of discrete loading 
is essential to understanding operation with a finite number of phases. 

It is shown that such generators may be modeled by an equivalent circuit 
containing a current source in parallel with a characteristic internal 
capacitance and a load impedance. The values of these circuit elements 
are determined by the excitation charge wave and geometrical factors such 
eis; load-to-charge wave spacing, charge-to-ground potential distance, and 
interelectrode spacing and number of phases in the discrete loading case. 

For the discretely loaded generator, the output theory is derived by assum¬ 
ing the form of the potential distribution along the loading electrodes. 

The discreteness requires that Fourier techniques be used and that the dis¬ 
tribution be treated as an infinite sum of spatial harmonics of the funda¬ 
mental wavenumber of the exciting charge wave. A six-phase, discretely 
loaded generator with an"electrode-to-interelectrode width ratio equal to 
1 . 5 ^ is analyzed in detail. It is shown that this generator can deliver 75^ 
of the power available from a continuously loaded generator having the same 
excitation and physical parameters. Data are presented to predict output 
power for discrete loadings, with the number of phases varying from five to 
42. A six-phase laboratory generator is used to determine the equivalent 
circuit element values. 

The theory predicts the cvirrent source mgignitude to within 6^. With the 
stray capacitance due to electrode structure quantified and included, the 
characteristic internal capacitance is predicted to within 5^ by the theory. 


Thesis Supervisor: James R. Melcher 

Title: Associate Professor of Electrical Engineering 
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CHAPTER I 


INTRODUCTION AMD BACKGROUND 

The past several decades have seen the development of electromechanical 
energy conversion machinery which has provided our industrial societies with 
a means for rapid technological growth. The scientists who experimented with 
electrical phenomena a century ago would not have predicted the posture of 
electricity in today's world, even in their most optimistic dreams. Since 
its implementation, electrical power has become an integral part of man's 
daily activity. 

Tne dependence of modem society on electric machinery is due to strong 
economic factors. Kinetic energy from a natural water fall may be converted 
cheaply and efficiently to electrical energy and then transmitted, with small 
loss, over great distances to industrial and population centers so as to pro¬ 
vide useful work. 

With relatively few exceptions, magnetic fields serve as the energy 

coupling mechanism for modem electromechanical energy conversion equipment. 

Until recently, the development of electric-field-based energy conversion 

machinery as a practical source of significant amomts of electrical power 

has received only passing interest. Historically, electric-field energy 

conversion devices were of great interest to those most intently involved 

in the development of electrical science during the nineteenth century. 

These machines were for the most part direct current generators with an 

energy storage capability, and were employed principally as experimental 

(9) 

laboratory apparatus. The twentieth century thus far has seen the Van 
de Graaff high-voltage direct current generator as probably the most signi¬ 
ficant development in electric-field energy conversion machinery. 
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The lack of development of electric-field machinery for power produc¬ 
tion has been due largely to severe limitations imposed by electric field 
breakdown strength at atmospheric pressure. Magnetic field systems are not 
limited by a breakdown phenomenon, and by employing them as the coupling 
medium, energy densities can be attained which are much greater than those 
attainable with electric fields before breakdown. 

With the magnitudes of the electric and magnetic field energy densities 

as a direct measure of their capacity to deliver power, a comparison of the 

(3) 

two types of machinery as energy converters may be readily accomplished 

1 3 

Magnetic field energy density: = flT 

Electric field energy density: Eg = — e Joules/m. 

Now, to compare the two energy densities, consider 

n = Eg/Eg 

To evaluate this ratio it is necessary to assign magnitudes to the mag¬ 
netic and electric field intensities. For field strengths normally attainable 
without the use of coolants or vacuum, systems, we have: . • , 

B =1 weber/m^ 

E = 3 X 10 ® voits/m 

Using these values and recalling that for atmospheric air 

y = UlT X 10 ~^ weber/amp-turn m 
e = 8 . 85 U X 10 ~^^ farad/m 






the relative energy density ratio may be evaluated: 
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n ^ 10** 

This calculation has exposed the advantage which has propelled magnetic 
field equipment into its present monopolistic position among low-frequency 
electromechanical energy conversion machinery. 

However, modem technology has made available high vacuTom systems and 
techniques for efficiently employing them which promote a questioning atti¬ 
tude toward the efficacy of the above energy density comparison. In a high- 

( 7 ) 

vacuum environment, it has been demonstrated that an improvement by a 
factor of 30 in electric field breakdown strength is not out of reach. 

In light of this information it is important to ask the question, **If 
given a very large electric field energy density, can it be practically em¬ 
ployed in a generator as the coupling medium for power production?" It is 
with the objective of answering .this question, at least partially, that this 
paper is being prepared.. 

The impetus for the research was provided by the recent work done inves¬ 
tigating electrofluid-dynamic (EFD) power. sources, conducted by the Continuimi 
Electromechanics Group here at M.I.T. In txirn, the need for such equipment 
has been supplied by the space industry and their applications for light¬ 
weight electric power generators aboard space vehicles. Electric-field-based 
generators offer a distinct advantage in power delivered per unit weight over 
their magnetic field counterparts. 

An EFD generator, characterized by a high voltage, low current output, was 
proposed and constructed as the subject of two recent M.I.T. Theses^^^ 

The scheme employed a high velocity air flow with entrained particles which 
were sinusoidally charged in a corona exciter. The flow, a traveling wave 
of charge, passed from the exciter through a long, slender channel surrounded 
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by a conducting medium. The charge wave interacted synchronously with the 
conductor, converting kinetic energy from the particle flow into electric 
energy which was dissipated as heat. The energy conversion interaction used 
in this generator is the electric field analogue of the interaction which 
converts kinetic energy to electric energy for magnetic field synchronous 
machines. 

The synchronous electric field generator proposed and built during the 
research for this paper resisted from an attempt to model the charged par¬ 
ticle flow realized in the EFD generator in order that loading character¬ 
istics could be studied. Details of the generator built aj:e set forth in 
Appendix B. A standing wave of charge density on a rotating disk was used 
to approximate a traveling charge wave passing through a long slender■chan¬ 
nel. Synchronism was obtained by establishing an integral number of wave¬ 
lengths of charge distribution.around the disk periphery. 

Other attempts to build electric field synchronous generators are not 

( 7 ) 

unknown. However, these have been voltage-constrained machines depending 
on a spatial rate of change of capacitance to provide the mechanism for 
energy conversion. These types of synchronous machines eire the electric field 
analogue of the synchronous reluctance machinery in the magnetic field system. 
The EFD generator discussed and the one built for this thesis research are 
charge-constrained synchronous generators requiring a spatial gradient in 
the exciting charge density, and no spatial capacitance gradient, for energy 
conversion. The magnetic field analogue for this machine is what is known 
simply as a 'smooth-air-gap synchronous generator'. 

It is interesting to note some of the salient advantages that the charge- 
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constrained machine possesses over the voltage-constrained machine. 

The electrode hlade structure in the voltege-constrained machine under¬ 
goes very rapidly alternating electric pressure forces which present 
designers with the possibility of an instability occurring in one of sev¬ 
eral possible modes of vibration. The severity of this problem is obvious 
when the blade spacing is of the order of a millimeter and the rotational 
speeds approach 30,000 RPM. The torque-time history for such a machine 
shows rapid fluctuations due to the time-varying electric traction as the 
blades move in and out of each other's influence. Both of these problems 
impose severe limitations on the feasibility of high rotational velocity 
operation, which is essential to the production of large amounts of power. 

In contrast, the charge-constrained synchronous generator does not 
require a spatial gradient in capacitance, hence removes these limitations 
on the feasibility of high speed operation. For this machine, the elec¬ 
trical traction force on each part of the rotor (and therefore the mechani¬ 
cal torque) is constant in time. Also, the severity of the alternating 
pressure forces acting on the electrode structure is greatly diminished 
with several electrodes per wavelength of induced voltage. For continuous 
loading, as is discussed in Chapter II, this effect is even further reduced. 

In summary, this research concerns itself with a charge-constrained 
synchronous generator placing particular emphasis on providing a theory by 
which the output of such a generator might be predicted from its physical 
parameters and charge wave description. An actual generator was built to 
allow for a direct comparison of theory and experiment. Through this com¬ 
parison, it is hoped that a step forward in the development of a new scheme 
for the generation of electrical power has been made. 
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CHAPTER II 

CONTINUOUS LOADIMG 

Introduction; 

A detailed analysis of the continuously loaded, charge-constrained 
synchronous generator can provide an excellent insight into the energy 
conversion interaction involved. This development will serve as the basis 
for subsequent analysis in Chapter III of the discretely loaded generator. 
The geometry of the continuously loaded generator will be modeled as shown 
in Figure 2.1. The excitation source for this generator is a sheet of 
sinusoidally distributed surface charge density traveling along the chan¬ 
nel centerline at a fixed velocity. To an observer fixed in the frame of 
the load, this would appear as a traveling wave of charge. The power out¬ 
put of the generator will be measureable in terms of resistive heating in 
the conducting load. In this chapter, the source and load will be reduced 
to lumped circviit elements normalized to channel load interfacial area. 

This equivalent circuit approach will prove to be a powerful technique for 
studying such a generator and will be extended to the discretely loaded 
case in Chapter III. 

Governing Equations: 

From Fig. 2.1,this problem involves two distinct regions, the load and 
the channel. The most direct method for determination of the fields is to 
assvime a potential distribution at the load-chsmnel interface, solve the 
bulk equations independently in each region, and then match solutions at 
the interface. Because of the symmetry of the problem about the y=0 axis, 
it is apparent that the problem need be solved only for the upper region, 
with solutions for the lower region taking the same form. 












xU 



Figure 2.1 


Continuously Loaded Generator 
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From Maxwell's equation, the governing equations for the two regions 
of interest are: 


(a) 

Region 

(1) 

0 

0 _< y _< d 

(b) 

Region 

(2) 


d. <y < d + A 


Since the potentials in both regions are established due to the influ¬ 
ence of the traveling wave of charge, they will be assumed to be traveling 
waves in form also, with the same wavenumber as the charge wave. Potential 
variations in the z direction may be neglected, as the channel is assumed 
slender and wide such that 2d/W« 1. On this basis, the potential, distri¬ 
butions must take the form: 

♦ = Re My) eJ<"‘ ' (2.2) 

Substitution of the assumed form of the potential solution (Eq. 2.2) into 
the governing equations (2.1) yields for each region: 

- k^i(y) = 0 (2.3) 

The potential solutions therefore are: 

A* 

(a) cosh ky + sinh ky 

A 

(^) <l>j^(y) = cosh ky + sinh ky 

Boundary Conditions; 

Evaluation of the four constants in equation (2.1*) will complete the 
determination of the potential distribution. This requires four independent 


0 ^ y£ d 

(2.1t) 

d < y < d+ A 
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boundary conditions on the fields. By specifying an assumed potential dis- 
P^lbuticn along the interface as Indicated on Fig. 2.1. and by specifying 
the location of tero potential at the outer boundary of the conductor, 
three boundary conditions are available: 


(a) 

(b) 

(c) 


= V 

c|)j^(y=d + A) = 0 


(2.5) 


The last boundary condition is obtained from the charge sheet. At the sheet, 
the normal component of the electric field is e,ual to one-half the magnitude 
of the charge density divided by the dielectric constant. The factor of one- 
half is due to the symmetry of the channel; half of the charge excites the 
upper region and half excites the lover region. This boundary condition 

becomes: 

(^) 2E, 



Q 


Solutions: 

Applying the boundary conditions of equation (2.5) to potential solu 
tions expressed in equation (B.U) leads to the follo.ing matrix equation for 

the unknown constants. 


cosh kd sinh kd 0 0 




A. 

V 

A. 

0 1 0 0 




- '^/2e k 
o 

0 0 cosh kA sinh kA 

0 0 1 0 




0 

A. 

V 


( 2 . 6 ) 
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Evaluation of the constants and substitution into Equation (2.4) com- 
plately specifies the fields in terms of the charge density and the assvimed 
potential distribution on the interface. 


(a) <|)Q(y) 


■ fe 




cosh kd 2e k 
o 


tanh kd 


cosh ky - 


A 


2e k 
o 


sinh ky 


(2.7) 


(b) ^ cosh k(y-d) - V coth kA sinh k (y-d) 


Equivalent Circuit; 

Using equation (2.7) it is now possible to develop an equivalent cir- * 
cuit which characterizes the channel and load interaction. First the load 
is characterized by calculating the total current density, conduction current 
plus displacement current, which flows into the load region normal to the 
interface as a function of the assumed potential distribution. This load 
current density is a traveling wave with the same wavenximber as that of the 
potential and its complex amplitude is given by: 

3 D 


(a) 


u 

/y=c 


= aE 


yy 5 1 

y=d / y=d 


( 2 . 8 ) 


In terms of the potential distribution of equation (2.7 b), 


(b) 


/ y=d ^ y=d 


Performing this calculation and defining equivalent circuit elements completely 
characterizes the load region. 


Jt = [jwC^ + G_]V amps/m^ y=d 


where; 


(a) 

(b) 


^L = 


ke^coth kA farads/m^ 

ii 

ak coth kA mhos/m^ 


(2.9) 

( 2 . 10 ) 
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The channel may now he characterized in a manner similar to that used 
for the load to obtain the current density normal to the interface out of 
the channel region. 


^ A. 


Jq = [- + Jg] amps/m^ y=d 


( 2 . 11 ) 


where an equivalent capacitance and source have been defined. 


(a) C_ = ke tanh kd farads/m^ 
C o 


( 2 . 12 ) 


(b) J, 


.1 cuQ 


2 cosh kd 


amps/m^ 


The equivalent circuit may now be obtained by recognizing that there is no 
time rate of free surface charge at the interface, and therefore the current 
densities out of and into the interface,as represented by equations (2.9) 
and (2.11) respectively,must be equal. 

^ /N ^ 

[ + G^]V = - jo) CgV + Jg (2.13) 

The equivalent circuit for the generator which represents the load-channel 
interface as a terminal pair is sho™ below in Figure 2.2. 



Fi gure 2.2 
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Pover Output Calculation; 

The computation of the time average pover delivered to the load or. a 
per square meter of chaJinel-load interfacial area basis is straightforward, 
making use of the equivalent circuit of Figure (2.2). The power output 
found in this way is the power delivered to the upper load region only* it is 
multiplied by two to account for the lower load. The total power output is: 


P> = Re 


r j j* 

L J 


watts/m^ 


( 2 . 11 +) 


In terms of known quantities: 


J = 


__ 

ja)[C^+ C^] + 


amps/m^ 


(2.15) 


therefore: 


< P > = 




watts/m 


(2.16) 


Power Optimization: 

The power output expression of equation (2.l6) can be plotted as a func¬ 
tion of load conductance for fixed channel and load geometry and fixed charge 
flow condition. 



- G 


L “ 


Fig. 2.3 
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From Fig. 2.2, it is obvious that there exists some optimum conduc¬ 
tance for which the maximum amoiint of power is delivered across the 

interface to the load. Tc locate this optimum, maximize Equation (2.l6) 
over G^t 


9<P> 


= 0 


(2.17) 


Ceurrying out the calculation specified by Equation (2.17) yields: 


’^OPT ~ mhos/m^ 


( 2 . 18 ) 


The optimum power is obtained by substitution of Eq. (2.l8) into (2.l6): 


<P> 


u 


M 


OPT 


8[C„ + ]cosh2 kd 

L/ Ij 


watts/m^ 


(2.19) 


It is important to note that Eq. (2.19) is only a limited optimization. 
In general. 


<P> = f ( Q,to,a,d,A,£ ,k) 

O Ij 


( 2 . 20 ) 


This function may be considered as an eight-dimensional vector which 
traces some complex surface as the values of the independent variables 
are ranged. The absolute optimum power can be found only be considering 
all eight variables. The preceding optimization ranged but one variable — 
conductivity. 

Another important factor to take into account when optimizing is the 
electric field breakdown strength.It will determine whether or not an 
optimum point is indeed physically realizable. 
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CHAPTER III 

DISCRETE LOADING 

Introduction; 

The continuously loaded, charge-constrained synchronous generator 
analyzed in Chapter II provided a simple mathematical model which placed 
particular emphasis on the seQ-ient features of generator operation. In 
terms of this model, the fundamental energy conversion interaction may be 
readily understood. However, as a practical matter, a generator of this 
type is of very limited value for use as a source of electrical power. 

Most equipment requiring electrical power for operation has discrete elec¬ 
trical terminals, and cannot be considered as a continuous type of load. 
Therefore, of more direct interest in determining the feasibility for power 
production would be an analysis of a discretely loaded generator. Rather 
than an infinite continuous loading, there would be discrete electrodes of 
finite width used to tap off the electrical power, which could then be 
transmitted to the locations where electrical machinery and/or equipment 
are employed. The number of phases would be determined by the nvmiber of 
electrodes per wavelength of exciting charge wave. 

An analysis of such a generator will be presented in this chapter, 
using the techniques developed in Chapter II for obtaining an equivalent 
circuit. The physical characteristics of the generator under discussion 
are as shown in Figure 3.1. 

Voltage Distribution Model; 

As in the case of a continuvan load, it is necessaiy to solve a two- 
region problem governed by Laplace's equation. The principal difference 
between the continuum and discrete loading cases is in the form of the 












Load 

Impedance, 



J (a)t + k x) 


Load 
Electrode 
y = - d 





• • • 


Q= Fe 

U = to/k. 


$ = Re / <1> e*' n 

n 


n = -«> 


y =,-(d + A> 




uo 

P 

I 


depth w 


Figure 3>1 


Discretely Loaded Generator 
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potential distribution along the electrodes on the channel boundary. Due 
to the finite electrode structure, it is incorrect to assume a potential 
vave with the same wavenumber as the exciting charge wave. Rather, spatial 
harmonics of the fundamental wavenumber are introduced and the potential 
distribution takes the form of a sum of an infinite series of traveling 
waves. This 8l11ows the problem to be analyzed by means of Fourier tech¬ 
niques. 

In order to proceed with a Fourier analysis, it is necessary to obtain 
a model for the potential distribution along the channel electrodes. The 
first and most obvious constraint on this distribution model is that the 
potential across each electrode must be a constant, since they are assumed 
to be perfect conductors. To determine the potential distribution between 
adjacent electrodes, to be used in the model, two possible courses of action 
present themselves. The first is to specify that there is zero net current 
density flow normal to the channel axis into the load channel interface. 

That is to say, there is no time rate of change of free charge on the inter¬ 
face in the interelectrode regions. The second course of action is to assimie 
the form of the potential distribution between electrodes, based on physical 
reasoning. 

With the electrode spacing assumed to be small compared to a wavelength 
of the exciting charge wave, a linear potential distribution between elec¬ 
trodes is felt to be a good approximation, and therefore the second course 
of action will be followed in the ensuing analysis. It is acknowledged that 
the first alternative would provide a more exact solution; however, the 
increased complexity of the problem would serve only to confuse, rather than 
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to expose the hasic principles involved in discrete loading. It vlll be 
shown in Chapter IV that the results predicted by this model are reasonably 
close to those realized from an actual generator. It should be noted here 
that the slope of the assumed linear distribution must be kept finite, as an 
infinite slope such as would resvilt from a stepped potential distribution 
* implies an infinite capacitance. Fortunately, as will become clear in the 
analysis, the capacitance is not a strong function of slope. 

Now that the general form of the potential distribution model has been 
determined, the pieces must be put together in such a way as to maintain the 
identity of each electrode in the solution. This is desirable, since, when 
using the generator, each electrode's output must be determinable in terms 
of its own loading and that of each of the other electrodes. This may be 
accomplished by isolating all electrodes of the same phase as a periodic 
potential distribution in space, which can then be Fourier-analyzed inde¬ 
pendently of the electrodes contained in the other phases. The total poten¬ 
tial distribution in space can then be obtained by summing these Fourier 
series over the number of phases per fundamental wavelength. The identity 
of each phase is maintained by assigning to it an index. The time frequency 
of each Fourier component in this analysis will be identical to the time 
frequency of the charge distribution wave. 

Consider the case where there are s electrodes per fundamental wave¬ 
length, then the potential distribution associated with the r'th electrode 
phase will be assumed to be of the following foiro: 
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c = a + t 
X = sc 


Figure 3»2 


V (x) = 

r 






(r-l)c ^ X < (r-l)c + 'b 


(r-l)c+b < X ^ rc 


re ^ X < re + b 


(3.1) 


A wavelength of potential distribution for the case of six (6) electrodes 
per wavelength would then have the following form. 
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Fourier Analysis; 

In terms of the preceding discussion, the potential distribution on 
the electz'Odes may he expressed as 




(3.2) 


The complex amplitudes are determined from the infinite sum. 


00 


, , 1 > - 
= 77 ^ ''r ® 

n= -oo n 


vhere 


k = 


2m 


n sc 


(3.3) 


(3.4) 


The Fourier coefficients for the series (3.3) are obtained from the follow¬ 
ing sum of integrals using Equation (3.1). 


(r--l)c +b 
f 

V 


V = 

r 

n 


rc 

r 


[x-(r-l)c]e dx + 


-JV 

V e dx 

r 


(r-l)c 
rc +b 
f 

V 


(r-l)c+ b 


(3.5) 


-Jk X 

r r-. 1 n ^ 

” Lb " X + rc J e dx 


r c 


Performing the above-indicated integrations yields the following Fourier 
coefficients for the r’th electrode: 


V r k c k 

-jr 


b 


I 2 . n 

4 e sin sin - 


f] 


(3.6) 













-18- 


Using this result together with Equations (3*2) and (3.3) completely speci¬ 
fies the complex amplitude of the potential distribution along the electrodes 
in terms of the assumed model and its phase amplitudes. 


v(x) 


s ^ 

E E 

r=l n= 


V 


scb 


k2 


f J^n^ 

l+e-T" 


Sin 


2 -2 
a sin J 


2 I J^n* 


(3.7) 


This completes the Fourier analysis and provides the basis for characterizing 
the generator by an equivalent circuit. 

Governing Equations; 

From Fig. 3.1 the symmetry of the generator channel about its center- 
line is apparent. It is therefore only necessary to solve the fields prob- 
■ lem in the upper channel and load region, since the lower solution must take 
the same form. The governing equations for the two regions are the same, as 
in the case of continuous loading. 


(a) 

<3 

II 

O 

0 ^ y d 

(b) 

= 0 

<il.y<_d + A 


In the preceding section the potential distribution along the load electrodes 
was analyzed as an infinite sum of harmonically related traveling waves. 

Since the potential in both the channel and load regions must equal the 
assumed distribution at the interface, it is efficient to think of the po¬ 
tential distribution in these regions as also being made up of an infinite 

sum pf traveling waves, harmonically related to the exciting charge wave. 

Each of these harmonics will therefore be required to satisfy equation (3.8) 

independently. That is to say, for both regions of interest the potential 

:rave will be of the form: 
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VO 

J(x ,y ,t) = Re Et (x,y,t) 


ns: - oo 


where 


J^(x,y,t) = $^{y) e n 


(3.9) 


( 3 . 10 ) 


For such an analysis, the governing equation (3.8) is also rewritten as: 


(a) 

= 0 

0 ^ y c 

(3.11) 

(b) 

o 

II 

> 

d ^ y <_ d + A 



Substituting the assumed form for the solution, equation (3.10) into the 
governing equation and solving determines the potential distribution in 
terms of four (U) vinknown constants. These constants must be determined 
for each harmonic. 


(a) 

C 


(b) 

L 



cosh kny 


sinh k^y 


cosh k V 

+ 

sinh k y 

“l 

n*' 


n*' 


( 3 . 12 ) 


Boundary Conditions: 


As with continuous loading,specification of the potential distribution 
along the interface and on the outside bovindary of the load region provides 
three (3) of the required four boundary conditions. These can be expressed 
for each harmonic using Equations (3.7) and (3.12). 


(a) 

(b) 


(y=<i) = V 
<IfaL(y=d) = 

As 

d +A) = 0 


(3.13) 


(c) 
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The foiirth bovindary condition is provided by the charge distribution, again 
as done in the continuous load case. However, only the fundamental harmonic 
is involved in matching the channel distribution to the exciting charge 
wave. 


(d) 


d'J'j^(y)' 

dy 


'y=0 




where: 



Solutions: 


Application of the hoiandary conditions expressed in Equation (3.13) 
to the governing equation distributions of equation (3.12) completely 
specifies the potential distribution in the load and channel regions in 
terms of the phase complex amplitudes. 



n 


Q 6 


-in 


cosh kjjd 


2e k 
o n 


Q6 


cosh k y 


*in 


2e k 
o n 


sinh k y 


(3.1U) 


- V [cosh k (y-d) - coth k A sinh k (y-d)] 
n n n n 


where; 


-J^n^c 


; .V V , r. 


jkna/2 ^ k c 


Sin sin 


m 2 10 n 


(3.15) 


(3.16) 


r=l 


as developed in the Fourier anedysis section. 

Equivalent Circiiit; 

The analysis involved to obtain an equivalent circuit for the discretely 
loaded generator is analogous to that followed in Chapter II for the continuum 
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loaded generator. To characterize the channel and load regions, solve for the 
net current density normal to the interface as a function of the potential 
distrib-ution along the electrodes. Then to include the effect of the lumped 
loading, as shown in Fig. 3.1j consider the 2- electrode and apply conser¬ 
vation of current. 


/V 



V Y 

z z 


Zc 



(3.IT) 


+ b n= - “ 


In the expression the current density is also represented as an infinite 
series of harmonically related traveling waves. The complex amplitude for 
each harmonic is determined by the lumped circuit characterization of the 
load and channel regions. 






+ C 


n 





(3.18) 


The source term is the same form as found in the continuous case. 

■^S ' 2 cosh 

The lumped capacitance characterizing the channel region is: 


(3.19) 


C = e k tanh k d farad/m^ (3.20) 

nc on n 

and that characterizing the load region: 

C ^ = e k coth k A farad/m^ (3.2l) 

nL o n n 

Interchanging the order of summation and integration in Equation (3•IT) and 
carrying out the integration yields: 
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g 


JuQ sin -jk (s,c- a/2) 


k cosh k d 

1 1 


r—1 "V ^ 


[C + C.][l-e 
nc nil-'*- 


n= - 00 


( 3 . 22 ) 


Equation (3.22) is the most general form of solution for the voltage on the 
electrode as a fiinction of its load admittance, and the voltages and 
load admittances of all other s — 1 electrodes. It is convenient to consider 
Equation (3.22) as a vector equation of the form. 






— A “ 

a a - a _ 


V 


D 

11 12 iS 


1 


1 



A 


A 

a 


V 


D 

21 

» 


2 

1 

= 

2 

1 

X 

1 


» 

A 


X 

A 

a a 


V 


D 

Si ss 


s 


s 

— - - 




— 


(3.23) 


or in vector notation: 


A. A. 

A V = D 


(3.2l*) 


where the elements of A are of the form: 


It 




f CO 

-f hr ^ ^ 

n=: -00 


8V 


schk^ 


-[tanh k^d + coth k^Aje 


jk^(il-r)c 


sin- 


k^a 


^ b 

sin sin 


k c 
n'' 


2 o 


(3.25) 


and of D 


jcoQ^sin k^ a/2 -Jk^(£c-^/ 2 ) 


k cosh k d 

1 1 


( 3 . 26 ) 
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Since the matrix A is square and therefore in general has an inverse, it 
is now possible to solve for each of the s electrode voltages and the 
problem is completely specified and solved. 

/V /V /V t 

V = A"‘ D (3.2T) 

This general case would he extremely difficult to solve and is of limited * 

interest. Of prime interest is the case of a uniform loading, all the 

electrode load admittances identical. In this case the magnitudes of all 

the V *s will be identical and they will be shifted in space from each other 
r 

by some phase angle. This can be shown to be, in fact, the case, by assum¬ 
ing electrode .voltages of this form, substituting them into Eq. (3.23) and 
noting that it then reduces to s identical equations. 

Assume: 

‘ 3 . 28 ) 

and Vjj = ^2* (3.29) 

Substituting this into Eq. (3.23) and noting that the doubly infinite sum 
on n is simply the sum of a complex number and its conjugate, the one re¬ 
maining equation for the voltage amplitudes is the following: 


juQ sin k^a/g 
kj cosh kjd 

y -Jkirc ~ 

o . n.._ J l6e [tanh k^d + coth k„A] , . , tr. . , ir^ 

— + JOdCo ^- ^2 -“— cos kj^rc sin kna/2 sin kjjb/2 sin k^c/2 

r=0 ® n=0 n 


(3.30) 


This represents the solution that the foregoing analysis had as its objective. 
It is interesting to note that, in taking the limit as b 0 and s -»• oo, this 
solution collapses to one identical in form to the one expressed by equations 
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(2.9) and (2.10). 

Equation (3.30) can be represented as an equivalent circuit which 
characterizes the discretely loaded generator in lumped elements on a per 
meter of channel depth basis. 



where: 


S”1 ^ 

^ X l6e [tanh k d+coth k A] , . k a . k t . k c 

= £ s > —^- ) _S_cos k^rc sin n sin n sin n 

T / . set Z. ,2 ^ — T 

^n 

farads 

juQg sin kia/2 

k, cosh k,d - <5.32) 


s-1 ^ “ 

Y l6e Y 

Z_J seb 
r=0 n=o 

and the source is defined: 


K = 


Pover Output: 

Using the equivalent circuit of Fig, 3.^, the output power for the dis¬ 
cretely loaded generator is readily calculated. Since only the upper half of 
the generator was ansilyzed, this power figure must he multiplied hy two in order 
to give the total power output. Consider pure resistive loading then: 
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P ^ = „2 . 27—2 watts/electrode 


Power Ootimization: 


( 3 . 33 ) 


As with the continaous load case, there is a matched loading which 
produces optimum power output. 

SpT ~ (3.3i») 

For this load, the optimum power becomes: 

IkI^ 

<P> = —— ' —watts/electrode (3.35) 

OPT SOJCrj, 

where K and have been defined in Equations (3.32) and (3.31). 

Comparison of Discrete and Continuously Loaded Generators: 


Having developed the theory for both the continuously loaded and the 
discretely loaded generators, it is now possible to make some comparisons 
between the two arrangements. In particular, it is of interest to examine 
the difference in characteristic internal capacitance and optimimi power 
outputs. In order to make such comparisons, it is helpful to rearrange 
Equation (3.30) into a form that is normalized to electrode area. 

, ' __ 

y _ _ 2 cosh k^d 

*^^i^*^cos k rcl tanh k d + coth k /Uf(n) 
n J L n n J 


'— + jwe^k 

ag o 


Z Z .n[e- 


n=0 r=0 


sin kja/2 


k.a/2 


-] 


( 3 . 36 ) 


where 


f(n) = 


. k a/2 
sin n 

k a/2 
n 


. k b/2 
sin n 

k b/2 
n 


. k c/2 
sin n 

k c/2 
n 


The sum over r in the denominator of this equation may be rewritten and the 
result used in defining a characteristic internal capacitance for the discretely 
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loaded generator. 

°° s-1 

= e k, / / 2n[cos k rc cos k rc][tanh k d + coth k A] f(n) farads/m 

D o i_j n n n 


n=0 r=0 


(3.37) 


Recalling from Chapter II that the characteristic internal capacitance for 
the continuously loaded generator is given by: 

= e k [tanh k d + coth k A] farads/m^ (3.38) 

Cox X 1 

a ratio which compares these capacitances is defined: 

CAPR = (3.39) 

Before making a comparison of optimum power outputs, note that the total 
channel area over which power is delivered to the loading is not the same 
in the two types of generator, and the difference must he taken into con¬ 
sideration- The amount of this difference is reflected in the ratio of 
discrete electrode width (a) to the fundamental wavelength divided hy the 
number of phases (c). 

a = a/c (3.Ho) 


Taking into account the relative area factor, a ratio comparing the optimum 


power outpute may be defined using the optimum power output expressions for 
the discretely and continuously loaded generator as given by Equations (3.33) 


and (2-l6) respectively. 


FO\m 


'^^^OPT DISCRETE 

CONTINUOUS 



(3.1+1) 
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Due to the double sums involved in the above ratios, the computer was 
used to make sample calculations for curve plotting, and the programs 
designed and employed for this purpose are contained in Appendix C. 

Curves 3.1 and 3.2 demonstrate how the ratios : expressed in Eq. (3.39) 
and Eq. (3.^1) vary according to the number of phases used per wavelength 
and the degree to which the discrete electrodes occupy the channel inter¬ 
facial area. It should be emphasized that these curves are only as good 
as the model from which they were obtained. For this reason, they should 
not be expected to give valid results at the extremes of the ranges. 

The range of a studied is felt to correspond to that range over which 
the discrete generator model is valid. Frcan curve 3.1, note that CAPR is 
a maximum overall s for a - .5, and decreases symmetrically for a greater 
or less than this value. 

Curve 3.2 demonstrates the effects of a and s on output power. It 
appears from these curves that large a and s would be the most desirable 
arrangement. The power ratio goes asymptotically to ct as s is increased 
to large values and is always less than ot due to the area difference 
between discrete and continuous loading. However, as a and s are increased, 
the spacing gets very small and electric field breakdown strength becomes 
of critical importance. 

Using curve 3.2, one is led to make the optimistic judgment that dis¬ 
cretely loaded charge-constrained generators compare very favorably with 
continuously loaded ones on the basis of power output. An energy conver¬ 
sion efficiency which quantifies the degree to which the maximum amount of 
power available is removed from a discretely loaded generator as a function 
of tt and s is displayed on curve 3.3. The information used to plot this 
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curve is contained in curve 3.2 and is not a reflection of mechanical, to 
electrical conversion efficiency. All generators considered are assumed 
to be 100 per cent efficient in this regard, as no losses have been included 
in the model. 

Experiments on an actual generator are now in order so that the theory 
on which the preceding calculations are based may be verified. 
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CHAPTER IV 

EXPERIMENTAL COMPARISOH WITH THEORETICAL PREDICTIONS 

Introduction; 

In Chapter III, a general theory was developed for determining the 
loadihg characteristics of the charge-constrained synchronous generator. 
Modifications to this theory, to account for the geometry of the experi¬ 
mental generator, have been developed and are presented in Appendix A. 

It is now possible to present a theory which applies strictly to the 
laboratory generator and to make predictions concerning its output 
characteristics. 

From Equation (3-30), the output voltage for any electrode in the 
uniform loading situation is: 


V = 
o 






scb 

r =0 n = 0 


I 


( 1 ^. 1 ) 


f(n) 


n 


where: 


-/ \ . . k a/2 sin ^n^/2 sin ^n“/2 (k,2) 

f(n) = cos knrc sin n ' ' / \ / 


In Appendix A, the modifications for D, C_ and C 
into account the character of the actual channel, 
term, using (A.11), was: 


are developed, taking 
L 

The result for the source 


£ = . 1 ( 0 $ 2 sin ^ 1^2 _ 

k^cosh k d[l + e^fj tanh k^ d] 


(lt.3) 


The modified capacitance for the channel is: 

e k [e f + tanh k d] 
p, _ on r 1 _ n 

Dq “ [l + e^f^ tanh k^^d] 


(a.10) 
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t 


I 


and for the load region, the capacitance becomes: 



e k [coth k 5 + e ] 
_E_E_2_I- 

[l + e coth k_<S] 


(A.5) 


where the following definitions have been used: 

1+e tanh k (h-6)tanh k^5 

f = -- 

e tanh k (h -5)+ tanh k 5 
r n n 


and 


(A.12) 


(A.6) 


For the laboratory generator as described in appendix B, the following 
geometrical factors apply: 


a) 


e 

r 


h 


b) 

kn6 . 

= .213 n 

c) 

k^(h 

-6) = 

3.87 n 

d) 

k d : 
n 

= .177 

n 

e) 

IT : 

= 33.8 

n 

n 



f) 

Area 

of electrode = 


Short Circuit Current Measurement and Charge Prediction: 

For the purposes of this chapter, consider the following circuit 
model for each generator electrode and its load. 


I 
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Where the circuit quantities are defined from Equation (U.l). and Appen¬ 


dix A. 


I = g D' 


and the internal capacitance. 

s -1 


n =6 


^n 


f(n) 


(It.6) 


Now, using Equation (It.5) it is possible to make measurements of the short 
circuit current and predict the magnitude of the charge density on the disk. 
This prediction can then be compared with the actual disk charge density 
as measured with the eleotrometer circuit described in Appendix B. 

From Eq. (It.l) (It.5) and (It.6), we have: 


I 

o 


*0 + 


and for the case vhere 

/s 


^o/a)C^ »1 


(It.T) 


(4.8) 


the output current will, to a very good approximation, be the short cir¬ 
cuit current. 

The inequality of Eq. (^.8) vras obtained 
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on the experimental generator by loading with a pure resistance that 

7 n 

Y = 2 X 10” ohms (^•9) 

o 

It will become clear later that this loading does indeed satisfy the 
condition imposed by Eq. (h.8). 

To determine the relationship between the short circuit current and 
the surface charge density, recall: 

1^1 _ 2g(jo sin ^ 1^/2 IqI _ (^.10) 

I ' ~ k cosh k d[l+e f tanh k d] 

1 1 r 1 


Since the electrode current is given by the integral of the current den¬ 
sity over the area of the electrode, and since the dimensions of the 
electrodes vary due to the cylindrical arrangement, it is appropriate to 
let 

Area of Electrode (U.ll) 

® ~ a 

where a is the electrode width over which the integration (3.17) was 
carried out. 

Using the numerical values presented in Eq. (4.4), a relationship 
for predicting the disk surface charge density from the source (i.e., 
short circuit) current is obtained. Evaluation yields: 


iQl 


3.4i|i| coul/m^ 


(4.12) 


The following short table compares the magnitude of the surface charge 
density as predicted from the short-circuit current and that measured 
with the electrometer to demonstrate the validity of equa- 
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tion (1+.12). Several of these comparisons were carried out over a 


period of time, and the results were consistent emd reproduceable. 



IqI predicted ( coul/m^) |q| measured (coul/m^ ) 


2U.0 X 10"® 


2U.8 X 10 “® 


15.0 X 10"® 


15.8 X 10"® 


23.9 X 10 -® 


24.8 X 10-® 


Table U.l 


A representative distribution of surface charge density as a function 
of angle around the source disk, as obtained by electrometer measure¬ 
ments, is given on Graph (U,l) 

Evaluation of Internal Capacitance Predicted by Theory : 

The theory predicts an internal capitance for the generator that is 
expressed by Equation (4.6). Since the terms of the infinite series gO 
as 1/n^ and the series contains terms of alternating signs, it is expected 
to converge very rapidly, making the evaluation of the internal capacitance 
through truncation efficient. In addition, further examination of the 
series over the r variable exposes the fortunate fact that terms in the 
double series have value only for particular values of n. To show this, 
interchange the order of summation in Equation (4.6) to get: 


00 



n=0 


(4.13) 
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where 





.i2ir (n+l )r -.12T T(n-l)r 
e_ s + e _ s 

2 


(U.lU) 


In Eq_. (U.lU) the exponential form of cos and the definition of 

Eq. (3.^) have been used. 

Examination of Eq. (U.lU) reveals: 


f (n) = 


- Ill'll ri“l n o o 

n=l; -= m; -= m m = 1,2,3 

s s 


otherwise (I of s roots of l) 


(U.15) 


This of course greatly facilitates the eveJ-uation of the internail capa- 
citguice predicted by the theory. For the case s = 6: 


C. 

1 




n = 1,5,7,11,13,17,19 - 


k c 
n 

2 


(l*.l6) 


Using the physical characteristics of the generator as presented in 
Appendix B and truncating the infinite series, one obtains for the inter¬ 
nal capacitance associated with each electrode, as predicted by the theory 


= U.7 X 10 farads 


Sources of Ci^pacitance Not Accounted for in Theory : 

The theory as modified in Appendix A does not fully account for all 
the capacitance that characterizes the generator. The most significant 
source of capacitance unaccounted for in the theory has to do with the 
finite dimension of the electrodes in the vertical coordinate and the 















-38- 


screws and nuts used to secure them to the load disk. 


Consider Fig. U.l. 



Figure U.l 


It would be extremely difficult to perform an exact calculation to deter¬ 
mine the amovint of this capacitance; however, approximate evaluations 
may be made readily. Consider the circuit demonstrating the effects 
about any electrode as in Fig. ^,2. 



Figure h.2 


Recalling the phase shift between electrode voltages, the added capaci¬ 
tance due to the structure is given by {h.lQ) considering the associated 
resistance ( 'Vl0^° ohms) as an open circuit. 

C = 2[C, (l - cos 5) + C, (1 - cos +-] 


(lt.l8) 
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Using the impedance bridge, it was possible to get an experimental 
measurement of the capacitance involved in Equation (4.l8) 

Cj = 1.5 pf 
= .75 pf 

C 3 = .25 pf 

C “ 0 

Using these results: 

C = 4.8 X 10 farads (4,19 ) 

s 


Unfortunately this capacitance is approximately equal in magnitude 
to that predicted by the theory as inherent to the generator. However, 
any stricture that might be devised to hold the electrodes would have 
some capacitance associated with it. With both the internal and struc¬ 
tural capacitances so small, the result of any attempt to further reduce 
them would be difficult to detect. The structure capacitance was inclu¬ 
ded in the equivalent generator circuit used to predict output. This 
total capacitance was considered as the characteristic capacitance of each 
electrode. 


C 


T 


C. + C 
1 s 

-12 

9.5 X 10 farads 


(4.20) 


Prediction of Matching Load 

The experimental generator contains five wavelengths of electrodes 
with six ( 6 ) electrodes per wavelength, as is described in Appendix B. 

To avoid the problems involved in loading each electrode as a separate 
source, the experiments were conducted with electrode outputs of the same 
phase connected in parallel. This required the attachraent of only six 
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loads to obtain characterizing output data. The modifications to the 
generator equivalent circuit are obvious: 

(a) 1-^51 

(b) 5 (J+.S!) 

✓N. ^ 

o o 

Using the generator so configured, that is, like phase assembled in paral¬ 
lel, it was possible to predict the loading for maximiim power output 
using the results of Chapter III and Equation (4.20). 

5 to mhos (4 • 22) 

17.9 X 10 ® mhos 
6 

56 X 10 ohms 

This represents the loading that is predicted by the theory for the 
generator to produce maximum power with five electrodes of the same phase 
wired in parallel. On a per-electrode basis, if each electrode was loaded 
individually, the required loading is as predicted in Equation (4.21): 

R = 280 X 10 ® ohms (^♦. 23 ) 

^OPT 

The effect of parallel operation in determining the matched load condi¬ 
tion increases the flexibility in application of such generators. 
Experimental Determination of a Matched Load: 

To obtain the output data from the experimental generator, for the 
purpose of identifying the actual internal capacitance, it was wired in 
parallel, as described in the previous section. 


(a) Y, 

^OPT 

and for 0)= 377 rad/sec 


(b) Y, 


OPT 


and 


(c) R, 


OPT 









- 1 + 1 - 


The scheme employed to determine the matching resistive load (and 
thus the total characteristic capacitance) was to examine a log-log plot 
of output current magnitude versus resistance loading magnitude. Con¬ 
sider the magnitude of the output current as a function of the source 
current, total capacitance and loading: 


111(1^ 


Hr + 

L CO c' 


Vi 


( 4 . 21 *) 


where k = 


taking the logj^ of both sides, 


1 

(OC 



constant - ^ log [R^ + K^] 
d 10 h 


(4.25) 


From this, two limiting cases become apparent as R. 
wide range. First: 


El» 


1 




-1 


is varied over a 


(4.26) 


and 



03 C 




0 


(4.27) 


By plotting log I ^ ordinate and log^^ as the ahcissa, the 

value of internal plus structural capacitance can he found graphically. 
This was done on graph 1+.2 As indicated, the characteristic capacitance 
found was: 


= 9*65 X 10 farads 


(1+.28) 


To verify this capacitance, several runs were made with the results con¬ 
sistent and reproduceable. A load capacitance was also added in parallel, 
with the characteristic capacitance and the results were almost exactly as 
predicted; see graph U.3 . 














The power output curve as a function of loading was- also calculated 
and plotted on graph U.h . Note that the maximum power output occurs at 
matched loading. 

Discussion of Results; 

The validity of the theoretical model assumed to characterize the 
charge-constrained generator has been demonstrated in the previous sec¬ 
tions. In particular, the magnitude of the equivalent current source may 
be accurately predicted by the theory for a known surface charge density. 

The experiments have shown the actual internal capacitance to be 
approximately twice that which would be predicted by the theory aJLone. 
However, this predicted capacitance did not include that capacitance 
associated with the structure of the generator. When this extra capaci¬ 
tance had been estimated and included a^ part of the characteristic in¬ 
ternal capacitance of the generator, then the theoretical prediction 
accurately matched the actual generator characteristics. However, it is 
felt that the most accurate method of characterizing a generator in terms 
of equivalent circuit elements, for application, would be to obtain a plot 
of equation (U.2U). The theory serves the purpose of providing an approxi¬ 
mation of the capacitance for design purposes, but due to the structure the 
actual capacitance would have to be determined after construction. This 
does not differ from the procedure followed for designing and characteri¬ 
zing magnetic field machines. 










Electrometer Reading - (Coulombs) 



Electrometer Reading 

Area of Probe 


coul/m^ 


Area of probe = 5.25 x 10*** 



Curve U.l; Disk Charge Distribution 


0 = Angular Position on Disk 
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Curve 1*.3: Determination of Characteristic Capacitance vith Capacitance in the 


Loading 


























. 1 + 6 ^ 


Power Output per 5 Electrodes Vs Load 


10 

+3 

H 

H 



I 

A 



Curve 4.U 


















CHAPTER V 


Conclusion ; 

The application of charge-constrained synchronous generators in 
the production of large amounts of electrical power is to a great extent 
dependent on the feasibility of power removal with a finite number of 
discrete electrodes. That is, such generators must have output terminals 
with a finite number of phases through which they can deliver power. 
Otherwise, they are strictly limited to continuvim type loading. 

In Chapter III a discretely loaded generator with an arbitrary nvim- ‘ 
ber of phases was analyzed. By modeling the potential distribution between 
adjacent electrodes as a linear fvinction, it was possible to obtain a ter¬ 
minal pair representation for each electrode's output. The equivalent 
circuit for each electrode was shown to be accurately represented by a 
current source in parallel with a characteristic internal capacitance. 

It was then demonstrated in Chapter IV that the magnitudes of both of 
these equivalent lumped elements was accurately predicted by the theore- 
ticeil model. The need for a means of obtaining an accurate estimation of 
that part of the internal capacitance due to the physical structure of 
the load electrodes was demonstrated in the experimental section of this 
research. 

Of prime interest in the use of discretely loaded generators is the 
amount of power diminution that occurs comparative to the case of contin¬ 
uous loading. It has been shown that this penalty need not be severe if 
the electrode spacing and number of phases are carefully chosen. With the 
power output per fundamental wavelength of channel length directly proper- 
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tional to the electrode area it hecomes important to make the inter¬ 
electrode spacing as small as field breakdown strength will allow. 

It is this factor which has the most significant effect in reducing 
power output when using discrete vice continuous loading. 

It has been shown theoretically that a discretely loaded generator 
can perform efficiently as an energy converter. However, in this regard 
there are several unanswered questions which may serve as the basis for 
further investigations in the development of a practical system. The 
ability of such a generator to run self-excited, from start-up, is of 
obvious interest. Mechanical energy storage schemes have been proposei^^ 
to accomplish just this, but is is yet to be demonstrated on an actual 
device. Questionable is the feasibility of constructing a generator with 
rotational speeds in excess of 30,000 rpm and providing it with a high 
vacuum environment so that high electric field intensities may be sup¬ 
ported. And lastly is the obvious problem as yet unconsidered, of 
adapting a six-phase plus generator to provide power for a three-phase 
world. Needed for this is a very efficient, lightweight phase transfor¬ 
mer or system of capacitors which would combine the phases. 
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APPENDIX A 


Modification to Theory for Experimental Generator 


Introduction 


For the generator built to provide experimental verification 
of the theoretical model considered in Chapter III, it was necessary to 
make some changes to the load and channel arrangements. Practical con¬ 
sideration dictated that the channel be asymmetrically loaded and that 
plexiglas disks be used to support the electrodes and to carry the 
exciting charge distribution. These changes do affect the field dis¬ 
tribution, and therefore must be reflected in any theory applied to 
predict output for this particular generator. Conceptually there are 
no changes; however, the equivalent source and internal capacitance- 
lumped elements predicted in Chapter III must be modified. A cross- 
sectional view of the laboratory generator is given in Figure A.l, 
and demonstrates these changes. 


0) 

m 



Figure A.l 
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Since the voltage distribution along the electrodes is assumed known 
it is efficient to consider the modifications to the theory in two steps. 
First, consider the effect of the plexiglas layer above the electrodes on 
the characteristic capacitance of the load , as defined in Chapter III 
Secondly, consider the effect of the plexiglas layer below the electrodes 
and the asymmetric loading on the source term J^, and characteristic 
channel capacitance Cj^^^both as defined in Chapter III. These modified 
equivalent lumped elements should be used in making theoretical predic¬ 
tions for the experimental generator. 

Effect of Plexiglas Layer Ab ove Electrod e s 

The theoretical developments for this and the next section are car¬ 
ried out using techniques identical to those used in Chapter III. Any 
question concerning procedure should be referred to this chapter, as it 
is not felt to be necessary to go into the same amoimt of detail in the 
expositions of this appendix. 

The geometry of the problem being considered is shown in Figure A.2. 



Figure A.2 
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This is a two-region problem governed by Laplace's equation and there¬ 
fore the harmonic traveling wave solution will be of the form: 


(a) 'j' , = A , sinh k (y-d)+ B cosh k (y-d) d < y < d + 6 

n3 n3 n n3 n 

(b) $ = A e 

ntf n 4 

where the boundary condition of zero potential at y = oo has been used 
and the solution axes have been shifted to inake calculation more efficient. 

The three remaining boundary conditions which specify this problem 
are given by 


(a) 

II 

d+6) = 

/\ ■ 

$ (y = d +5) 
n 4 

(b) 

e V $ 

\ 

e V $ \ 


p y n3 

) 

0 y n**) 


/ 

' y=d+5 

/ y=d+ 6 


(c) $ (y = d) = V 

n 3 n 

Using the boundary conditions (A.2) to evaluate the unknown constants of 
(A.I) completely detemines the fields. Using these results, it is now 
possible to obtain the modified expression for the load capacitance 
by considering the displacement current density noimial to the channel 
axis at the electrode interface 


and 


Completing these 


J = -jcoe [V $ ] 

n P ' y 


C = - e V f 

nL p y n3 , 

y=d 


calculations yeilds the modified canacitance, 


(A.3) 


(A.U) 
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where 


C* 


'n. 



2 


(A.5) 


(A.6) 


Effect of Asymm etric Loading and Ple xi glas Charge Disk: 

Evaluation of the effects of asymmetric loading and the plexiglas 
charge disk requires the solution of a three-region fields proolem which 

is governed hy Laplace*s equation. The geometry of the problem is as 

shown in Figure A.3. 

The Fourier component solutions of Laplace’s equation in the three- 
regions of interest will take the form: 

(.a) $ = A cosh k (y -d) + B sinh k (y- d) 0 ^ y ^ d 

' ni ni n ni n 

(b) <41 = A cosh k y + B sinh ky -6<y<.0 (A.?) 

' ' n2 n 2 n*^ nz n’' — 

(c) 4> = A cosh k (y + 5) + B sinh k^(y+5) -h£y^-5 

n3 n3 n ns n 

where again the solution axes have been shifted to make calculations more 
efficient. 

The six (6) boundary conditions required to evaluate the unknown cons¬ 
tants of equation (A.7) are given as: 


(a) =<i) = v^ 



Q 


(c) (y=0) = (y=0) 


(d) $ (y= - 5) = (y = -5) 

Uz n 3 
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Asymmetric Channel vith Plexiglas Disk 


y= d 


y=0 


» _ J (wt - k_x) 

$„ = Re V e'' n ' 

n n 


Q 


Q = Re Q e 


j(aJt- k| x) 


; ,////<pV/////^/7/7777777 

9 = 0 

© 


y=-h 


<J> = 0 

D 


Figure A.3 
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(e) e V $ 


= e V 




(f) * (y =- h)= 0 

n3 


Note that the displacement current density normal to the channel axis 
at the electrode interface is given by: 


J = - jtjoe k B 

n ^ o n m 


(A.9) 


It is therefore only necessary to determine the unknown coefficient 
in equation (A.7*a) using the boundary conditions of equation (A.8). 

Determining and using it in equation (A.9) specifies the modified 
source and channel capacitance terms for the physical arrangement 
of interest. These modified terms are: 


and 


where: 


nC 


= 


^*1 = 


£ k [e f + tanh k d] 
o n r 1 _ n 

[l + e f, tanh k d] 
r 1 n 

Ja3Q 


cosh k d[l + e f tanh k d] 
n r 1 n 

1 + r tanh(h - 6)tanh k 5 
^ ____ n 

e tanh k (h-6)+ tanh k 5 
r n n 


farads/m^ 


amps/m' 


(A.10) 


(A.ll) 


(A.12) 


These results combined with results of Chapter III provide the theory for 
the actual generator built for experimental measurements. 
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Experimental Generator: 


General Generator Arrangement: 

An asymmetrically loaded, charge-constrained synchronous generator 
was designed and built as part of this thesis project so as to provide 
verification of theoretical results. A sketch of the device built is 
presented in figure Rl. 

The physical measurements of this generator are given below and should 
be used in any application of theory to predict output. Note that a mean 
wavelength, and therefore mean wavenumber, is used to characterize the 
load section. This is necessary since the wavelengths at the inside and 
outside tips of the load electrodes are different, due to the cylindrical 
geometry. Since the theory considers a rectangular channel this consti¬ 
tutes an approximation; however, the results as shown in Chapter IV indi¬ 
cate that it is a good one. 


X 

1 

=: 

18.6 X 10“^ 1 

k 

1 


33.8 m"l 

6 

= 

.6 X 10 ^ m 

d 

n 

.5 X lO"^ m 

(h -6) 

= 

10,9 X 10 ^ m 

A 

electrode 


11.5 X lO"” m 

a/c 

= 

.61 

0) 

e 

= 

60 cps 

0) 

m 

= 

10 cps 

e 

r 

=; 

U.O 

e 

o 

=: 

8.854 X 10 "^^: 


m 
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Plexiglas was chosen as the material to be used for the structural sup¬ 
port of the load electrodes because of its inherent attributes such as 
easy workability and transparency and its electrical, properties are known. 
The plexiglas was chosen for the charge disk material principally for its 
extremely long characteristic relaxation time for bulk free charge. 

Charges placed on the disk by rotating it through the corona exciter sec¬ 
tion remain in place long enough to easily complete a full set of experi¬ 
mental measurements. The charge disk can be fully energized in less than ^ 
10 sec. and then the corona section may be turned off to eliminate noise 
in the electrode output signals. 

In order that the corona exciting field would be in synchronism with 
the charge wave established on the disk, it was required that disk rota¬ 
tional frequency be divisible into the exciting voltage frequency an 
integer number of times. For the laboratory generator: 

^electrical _ 60 HZ _ ^ 

^mechanical 

This establishes six wavelengths of charge distribution around the peri¬ 
phery of the disk. The disk speed was obtained by using an l 800 rpm 
synchronous motor driving the shaft with a toothed, flexible belt and a 
set of 3:1 speed-reducing pulleys. With both the motor and corona being 
driven from the same electrical mains, the maintenance of synchronism was 
not a problem. Electrodes were arranged on the load disk at six electrodes 
per wavelength of charge distribution. Although there were six full wave¬ 
lengths of charge distribution about the load section circumference, only 







-57- 


Appendix B 

five wavelengths of load electrodes were installed. The sixth wave¬ 
length was left clear to accommodate the corona source. 

Measurement of Disk Charge Density; 

An important variable in describing the generator output is the 
magnitude of the charge distribution wave on the disk. It is not pos¬ 
sible to gat an accurate measure of this charge magnitude simply by 
knowing the magnitude of the corona voltage. A direct measurement of 
the distribution was obtained using a Kuthly Electrometer model 6lOB in 
conjunction with a flat plate sensing proble of known area. The probe 
was mounted in the location where the corona source is set during the 
disk charging cycle. Manually rotating the disk then puts the probe 
under the influence of the charge* the magnitude being a fimction of 
position. 

The electrometer is a dc instrument which measures the charge by 
integrating the current required to set up imaige charges for that charge 
influencing its probe. For this reason, it was necessary to make the 
rotation of the disk for charge distribution purposes extremely slow. 
Using this scheme for charge measurement gave consistent, reproduceable 
results which compared very well with theoretical predictions as shown 
in Table U.l. 

Pictures of the loading section with the electrodes loaded individ¬ 
ually Eind in parallel are given in Figure B.2. These pictures also show 
the electrometer probe and corona source as they appear when mounted. 
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h-H 

^2 


Electrode 


Rj = lU.O X 10“^ m 

R = 21.5 X 10"^ m 

2 

\ l > = 60"^ 

0 = - 6 . 1 ® 

-2 

a = 1.0 X 10 m 

1 

a = 2.15 X 10"^ m 

2 

g = 7.5 X 10~^ m 

cj electrical 


Charge Disk 


Figure B.l 
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Ca; 

View Showing Load Electrodes of the Same Phase Connected in Parallel. 
Note Corona Source in Position in Foreground. 


(b) 

View Showing Each Electrode with Individual Resistive Loading. Note Elec¬ 
trometer Charge Measuring Probe in Position in Foreground. 

Figure B.2 
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Computer Programs for Calcxilation; 

The programs C.l and C.2 can he used to calculate the capacitance 
ratio CAPR and power ratio POWR of Equations (3.39) and (3.^1), respec¬ 
tively. These ratios provide a comparison between discretely loaded and 
continuous loaded, charge-constrained synchronous generators. If it is 
desired to calculate both of these ratios at the same time, the programs 
may be efficiently combined to decrease the computer time required. 

The required inputs have, for the most part, been defined in Chapters 
I and II. 

KD = k d 

1 

KDELTA = k A 

i 

ALPHA = a/c 
s = number of phases 

The integer M sets the upper limit on the summation over n. This series 

converges rapidly and six contributing terras provide very good accuracy. 

Recall from Chapter IV that terms in this sum are contributed only for: 

, (n +l) j ^ (n -1) . ^ 

n = q. -- integer; - — = integer 

s s 

The outputs are: 

H = capacitance for continuous generator 

F = total vsdue of double sum over s and n 

K = CAPR or POWR (depending on program being used) 
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JPOWR MAD 

START PRINT COMMEIITS TYPE KD, KDELTA, ALPHA, S, M, GO$ 

INTEGER N,R,S,M j 

READ DATA 

W'R GO .L.0.5, TRANSFER TO END 
PI = 3.11*159 
KC = 2*PI/S 
G = 2.0/S 

H = TANH.(KD) + 1/TANH.(KDELTA) 

KA = ALPHA*KC 

KB = (1.0 - ALPHA)*KC 

F = 0.0 

THROUGH SUM, FOR N= 1.0 ,1.0 ,N.G.M 
A = H*(TANH.(N*KD) + 1/TANH.(N*KDELTA)) 

B = (SIN.,(N*KA/2 )*SIN. (N*KB/2)*SIN. (N*KC/2))/(N*N*N*KA*KB*KC/8) 
C = 0 

THROUGH RSUM, FOR R = 0, 1.0,R.GE.S 
D = C0S.(2*PI*N*R/S)*C0S.(2*PI*R/S) 

RSUM C = C+D 

SUM F = A*B*C^'G + F 

K = H*ALPHA*SIN.(KA/2)*SIN.(KA/2)*U/(F*KA*KA) 

P'T V, H,F,K 

VECTOR VALUES V = $3E12.6*$ 

T'O START 

END CONTINUE 

E'M 


Program C.2 








-62- 


Appendix 

JCAPR 

START 


RSUM 

SUM 

END 

R 1.033 + 


MAD 

PRINT COMi''IENT$ TYPE KD, KDELTA, ALPHA, S,M,GO$ 

INTEGER N,R,S,M 
READ DATA 

W'R GO .L.0.5, TRANSFER TO END 
PI = 3.li^l59 
KC = 2*PI/S 
G = 2.0/S 

H = TAHH.(KD) + l/RAlffl.(KDELTA) 

KA = ALPHA*KC 

KB = (1.0 - ALPHA)*KC 

F = 0.0 

THROUGH SUM, FOR N= 1.0 ,1.0,N.G.M 

A = n*(tanh.(n*kd) + i/tanh.(n*kdelta)) 

B = (SIN.(N*KA/2)*SIN.(N*KB/2)*SIN.(N*KC/2))/(N*N*N*KA*KB*KC/8) 
C = 0 

through RSUIvI, for R=0,1.0,r.GE.S 
D = C0S.(2*PI*N*R/S)*C0S.(2*PI*R/S) 

C = C+D 

F = A*B*C*G + F 

K = H/F 

P'T V, H,F,K 

VECTOR VALUES V = $3E12.6*$ 

T'O START 

CONTINUE 

E'M 

.283 


Program C.l 
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